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Nature exhibits not simply a higher degree but an altogether different level of com- 
plexity. The number of distinct scales of length of natural patterns is for all practical 
purposes infinite. 



O ■ Benoit B. Mandelbrot [11 
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Abstract 



This paper presents a comparative study of two families of curves in W 1 . 
The first ones comprise self-similar bounded fractals obtained by contractive 
^ ' processes, and have a non-integer Hausdorff dimension. The second ones are 

unbounded, locally rectifiable, locally smooth, obtained by expansive processes, 
and characterized by a fractional dimension defined by M. Mendes France. We 
present a way to relate the two types of curves and their respective non-integer 
dimensions. Thus, to one fractal bounded curve we associate, at first, a finite 
range of Mendes France dimensions, identifying the minimal and the maximal 
ones. Later, we show that this discrete spectrum can be made continuous, 
allowing it to be compared with some other multifractal spectra encountered 
in the literature. We discuss the corresponding physical interpretations. 

keywords: Mendes France and Hausdorff dimensions, self similar curves, expanded 
curves, multifractal spectra. 
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1 Introduction 

Fractal geometry focuses on the local non-smoothness of physical and mathematical 
objects or sets, by means of their fractal dimension, of which there are several formu- 
lations: box dimension, Hausdorff, etc. These dimensions represent the "amount of 
space" occupied by the set, measuring its degree of "wrinkledness" , when viewed at 
smaller and smaller scales. Within this fractal universe we will restrict ourselves to 
bounded self-similar curves F in M. n , of Hausdorff dimension, dimu(F), larger than 
unity. The fact that dimu{F) > 1 indicates that F has infinite length, and therefore, 
confined to its convex hull, it has to be infinitely folded or wrinkled, in a self-similar 
way. This characterizes the first family of curves which we call Tu- 

The second family of curves we deal with in this paper, Tmf-, are unbounded, and 
free to spread their infinite length all around the n-dimensional space. Each such 
curve T is locally rectifiable and locally smooth, and hence with dim#(T) = 1. M. 
Mendes France defined a fractional dimension dim^F for them [2]. This dimension 
"looks at" such a curve V from afar, further and further away, instead of "at smaller 
and smaller scales", as mentioned above. The idea is to "zoom out" and study the 
scaling properties of their lengths when growing to infinity. 

Curves F in Tu are obtained by means of contractive processes and contractive ratios, 
we will obtain curves T in Tmf by expansive ones. 

It is worth commenting that this approach is completely different from Strichartz's 
[3] reverse iterated function system which constructs a new limit fractal set with the 
same dimension as the original. Our approach is to ascribe, to each F, a gamut of 
unbounded curves, locally smooth, together with their corresponding Mendes France 
dimensions. 

The paper is organized as follows: In Sec. 2 we recall the definition of self-similar 
bounded curves F in Tu by means of their finite number of contractive ratios a^, 
% = 1 . . . , N; in Sec. 3 we do the same with the Mendes France dimension via the 
expansive ratios 1/af, in Sec. 4 we construct, for a given F, a curve r ai in Tmf 
for each aj, i.e. we ascribe to F a finite spectrum of dimMF(r ai ). We identify the 
maximal and minimal dimensions of this spectrum as the Hausdorff and the divider 
dimensions of F. We test the sensitivity of dim^F: a minute change in the value of a, 
implies a variation in the dimension dimjv/^(r ai ). In Sec. 5 we make continuous the 
discrete spectrum obtained in Sec. 4. In Sec. 6 we compare this continuous multidi- 
mensional or multifractal MF spectrum with three multifractal spectra encountered 
in the literature: a) the spectrum of Renyi generalized dimensions, b) the thermody- 
namical formalism (a,f(a)) with an appropriate measure, and c) the corresponding 
range [a min , a max ]. Sec. 7 summarizes the conclusions. 

2 Curves constructed by similarities 

Let A — A\, A 2 , ■ ■ ■ , A N+ i = B be N+l different points in M n , satisfying dist(A ; , A i+1 ) < 
dist(A, B) = 1, for all 2 = 1, . . . , N. Let $ : R n -> R n be N similarities such that 
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Si(AB) = AiAi+i, so S{ are contractions, and = dist(A i; A i+ i) < 1 their ratios of 
similarities or contractors, l^i^N. 

Let pi be the polygonal whose vertices are Ai, A 2 , . . . , Aw+i, so it is formed by the 
N segments AiAi+x, this is the first polygonal approximation of the curve F . We 
call pi the generatrix of F. The polygonal p 2 is obtained by replacing each segment 
AiA i+ i by its copy S,(pi), which has the same endpoints A4 and A i+ ±. So p 2 has 
N 2 segments, and p 2 = UiLi ^(Pi)) an d so forth. Assuming the polygonal p^-i has 
been constructed, we replace the segment AiA i+ i by Si(pk-i), obtaining a polygonal 
p k made of iV fc segments, such that p k = [jf =1 <Si(pjt-i)- It is proved in [4] that the 
sequence of prefractals {p k } converges, as k^oo, according to the Hausdorff distance, 
to a limit curve F, satisfying 

TV 



so F is invariant for the iterated function system (IFS) Si, . . . , Sn, it is infinitely 
wrinkled, and has self-similar structure. The well-known von Koch curve (where 
ai — 1/3, i = 1,2,3, 4), or the curve in Fig. 1, are examples in R 2 . 
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Figure 1: Self-similar curve F in M. 2 ; generatrix p\ with 01 = 02 = 03 = 1/7, 04 = 2/7 
and 05 = 4/7. 

Recall that an IFS Si, . . . , Sjy satisfies the open set condition (OSC) [5, 6] if there 
exists a non-empty bounded open set U C M. n such that 
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with this union disjoint. This criterion guarantees that the components Si(F) do not 
overlap "too much" . 

In what follows, we will call Tu the family of self-similar bounded fractal curves, F, 
constructed as described above, and satisfying the OSC. Under these hypotheses, it 
is known that the Hausdorff dimension of F is the unique value d— dim^(F) that 
satisfies the similarity equation [5] 

N 



i=i 

3 The Mendes France Dimension of Expanded 
Curves in W 1 

Let Tmf be the family of curves rcK" that are unbounded, locally rectifiable, 
and locally smooth, i.e., any arc of T has finite length. We will give an idea of the 
"fractional dimension" defined by Mendes France for this type of curves [2]. 
For a curve T&J-'mf, we fix an origin and consider the first portion Tl of V of length 
L. Let e > be given and let Tl(e) be the e-parallel body of Tl, also known as the 
e-Minkowski sausage of Tl 

T L (e) = [J B(x 7 e) = {xe W 1 / dist(x,r L ) < e} . 

Let Al be the diameter of the convex hull of T^. Then, the Mendes France dimension 
of a curve T is, by definition 

dim MF (r) = lim lim inf lo gf ^(e)) 

V ' e\0 L/oo log Ax V ' 

where ^(T^e)) denotes the n-dimensional volume of Fl(e). It can be proved that 
the limit, when it exists, has a value between 1 and n, and that it does not depend 
on e, so we can drop it and rewrite (1) as 

6xm MF [Y) = hmmf . (2) 

L/oc logA L 

This remark is very important, because, intuitively, it says that it does not matter how 
"fat" the e-Minkowski sausage is, but how the sausage "fills up" the space according 
to the development of when L grows. Therefore, we are dealing with a concept 
of dimension which does not look at the curve at small scales, as the Hausdorff or 
box-counting dimensions do; on the contrary, this dimension "zooms out", looking 
from afar at the behavior of the curve when its length tends to infinity. 

Notice that, when the limsup^^ exists and equals the liminfk/< 00 in Eq. (2), it 
suffices to consider a growing sequence such that A.L k+1 ^cAi t , for any constant 
c> 1; in particular Ai k =c k . 
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Figure 2: Construction of the expanded von Koch curve. 

4 The Discrete Spectrum of Mendes France 
Dimensions 

The two families, Tu and Tmf, have no curve in common since their curves have 
absolutely different geometric features; however, we will make a geometrical process 
that allows us to link curves of both families, and thereby to relate their respective 
dimensions. We briefly review the main concepts, geometrical ideas and theorems, 
given previously in a detailed form in [7] . 

To start, let us consider a strict self-similar FeJ-'h, where the generatrix p\ is made 
of by N segments of equal length a, < a < 1. For example, in the von Koch curve, 
N = 4 and a— 1/3. In the first iteration we construct p[, identical in shape to p\ but 
all segments having unit length: p[ is p%, expanded by a factor of l/a — 3. Next, p' 2 
is P2 expanded by l/a 2 = 3 2 , and so forth, as indicated in Fig. 2. Each p' k+l contains 
p' k : there is a process of inheritance which guarantees the existence of a limit curve 
T, continuous, locally rectifiable and locally smooth, unbounded, i.e. in ^Fmfi that is 
the "expanded" and "unwrinkled" version of F. Since all segments in p' k are of unit 
length, it is easy to see that fj, n (pk(e))^ieN k , and so 

dimMF(r) = u m !^!M£il = Um M£^) = J^H = dimH(F , . 

k-*oo logA fc fc^oo logl/a" —log a 

In the general case of some different values, a\ ^ 02 ^ • • • ^ ajv, as in the example of 
Fig. 1, we can enlarge by l/a±, . . . , 1/ajv expansive factors or expansors, and obtain 
up to N T's in J-mfi with different values of dim^F- Let us focus on Fig. 3, the 
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Figure 3: The first two steps (out of scale) of the construction of the limit curves 
r ai ,r a4 and r as as expanded versions of F in Fig. 1. Note that p^ inherits p^ (red) 
in each case. 

largest expansors are 1/ai, l/a 2 and 1/03. Take 1/ai for instance, and observe that 
in each polygonal p a k the shortest segment has unit length, and that p a k adds new 
segments, longer than those in p^}_i- so, in the limit curve T ai , all segments are longer 
than or equal to unity. If we take the smallest expansor 1/ajy (l/as in the same figure) 
we observe that in each p a k N the longest segment has unit length, and that p a k N adds 
new segments, shorter than the ones in p a k _ v So, the limit curve T aN has arbitrarily 
small segments. 

Finally, an intermediate expansor 1/a.;, i^l,N, like 1/04 in the same example, yields 
polygonals p a k with some segments smaller and others longer than those in p' k H _ 1 , so 
T ai has arbitrarily small and large segments. We have the following result [7] 

Theorem 4.1 Let F G Tb. with contractors < a\ ^ a 2 ^ • • • ^ ajv < 1/ the 

reciprocals 1/ ai ^ 1/ ' a 2 ~^- • • ■ ^ 1/ojv > 1 are the expansors constructing limit curves 
T ai , . . . , T ajv e Tmf respectively. Then 

dim MF (r ai ) < dim MF (r a2 ) < ■ • ■ ^ dim MF (r^) . 

We will call {diniMF(r ai ), . . . , diniMF(r ajv )} the discrete MF spectrum associated 
with F. 
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4.1 Identification of minimal and maximal diniMF 

We have also the following results [7] 

Theorem 4.2 Let F eJ^h, and T aN e Tmf the limit curve obtained by the smallest 
expansor Then 

dim MF (r^) = dim H {F) . 

Theorem 4.3 Let F eFh, and T ai G Tmf the limit curve obtained by the largest 
expansor 1/ai . Then 

dim MF (r ai ) = l+ V=\ J . (3) 

log — 

a 1 

We have identified the maximal dim^F of the discrete MF spectrum as the Hausdorff 
dimension of F. To identify the minimal one, let us recall the divider dimension 
or compass dimension of F: given e > 0, M e (F) is the maximal number of points 
Xq, Xi, . . . , x m in F (in that order) such that \xj—Xj_\ \ — e, for j — 1, . . . , m. If L e is 
the length of polygonal p £ that joins the Xj, then L e ^e (M £ (F)—l)p^eM e (F). Then, 
by definition of dim^,, we have [5] 

logM £ (F) t logL £ 
dim di „(F) = hm = hm H r . (4) 

e~*0 -log£ £^0.1 

to log - 

Next, take the polygonals pk which yield F and apply the e-dividing method, but to 
Pk (instead of F) with Ek = all in each step k (Ek \ when k — > oo, for a\ < 1). The 
length of pk is L fc , where L = X^ =1 ctj is the length of the generatrix p\. Since ai is 
the smallest of all ratios, then M £k (pk)~L k /a\ which implies, by Eq. (4), that 

lim fogMeM = Hm log(L fc M) _ 1 t logL 



e k ^o -\oge k k^oo -logaf log — 

ai 

which, together with Eq. (3), yield 

Proposition 4.1 Under the same hypotheses of previous theorems, we have diniMF(r ai 
identified with the concept o/dim^(F). 



4.2 Sensitivity of dimMF 

The following result [7] 

Theorem 4.4 Let F eTmf, and a\<ai, iy^l then 

dim MF (r ai ) < dim MF (r a o , 

implies that, should, e.g. a\ and, say, a-i be infinitely closer, still the Mendes France 
dimensions of T ai and T a2 would differ. In other words, T ai has only segments larger 
than unity and arbitrarily large, a2 > «i implies: some small and smaller segments 
will be introduced in T a2 . Should these "wrinkles" be arbitrarily small and difficult 
to "see", still they would increase the dim MF . 
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5 'Continuization' of the discrete MF spectrum 

Indeed, the discrete MF spectrum {dim M p(T ai ), . . . , dim MF (T aN )} can be made con- 
tinuous, and in order to demonstrate it, we start with a simple example: only 
two contractors, 0<6<a<lasin Fig. 4 (where 6 = 1/4 and a = 1/2). In 
this case, the spectrum has only two dimensions {dimj\/ j p(r fe ), diniMF(r a )}, with 
dim MF (r 6 ) <dim M F(r a ). T b is the very stretched version, so that it has all segments 
larger than or equal to unity; whereas r a , although expanded, is a more wrinkled 
version since it has segments arbitrarily small. 




a 

u 1 F 

Figure 4: Generatrix pi with 6=1/4, a — 1/2 and the yielding fractal curve F. 

Hitherto, one ratio is chosen and inverted to make the expansion process in each step. 
But there is no reason for not using the inverse of some other ratio, or all of them, in 
the same process, selecting one for each iteration. 

For example, choose the expansor 1/6 = 4 for the odd steps k, and 1/a — 2 for the 
even ones. We take p\ the generatrix of F in Fig. 4, and setting the origin at the left 
endpoint (only for clarity, not for necessity) we generate p\ four times longer than pi 
(Fig. 5.1). This is equivalent, in this case, to stretching the interval [0,1] to [0,4], 
plus adding a "square hump" keeping the shape of p\. Next, we take 1/a = 2 as 
expansor, enlarging p[ by 2, but setting the origin at the right endpoint of p[ (again 
for simplicity), and adding square humps in a proportional way, obtaining p' 2 (Fig. 
5.2). Polygonal p' 3 is got from p' 2 , 1/6 = 4 the expansor, left endpoint the origin, 
adding the humps... and so on (Fig. 5.3-4). There is inheritance: p' k+ i contains p' k , 
which guarantees the existence of a limit curve T. In this case, after k expansions, 
we have the same number of segments and square humps that we would have by 
expanding only by 1/6, and starting from the left, or by 1/a, starting from the right. 
But, expanding k times by 1/6 would yield a p'\ with a diameter A b k = (l/b) k = 4 k . If 
we had proceeded using 1/a, A k would be (l/a) k = 2 k . Polygonals p'\ and p'\ would 
be identical in shape, but with different scales; with the same number of segments 
and humps, but of different diameters, 2 k and 4 fc : i.e., they would be similar. Instead, 
if we expand k/2 times by 1/a and k/2 times by 1/6 no matter in which order, 
we will obtain p' k with identical number of segments and humps as p' a k and p' b k , but 
now with a diameter 

A k = (l/6)t(l/a)f = 4*2* = ((l/6)^(l/a)^) fc = {^/W>)\ 

a value strictly comprised between l/a k and l/6 fc . This results exactly as if we had 
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Figure 5: Polygonals p\, p' 2 , p' 3 and p' A of expansion process corresponding to Fig. 4- 
Each new polygonal (in red) inherits the preceding one (in black) (steps 3) and 4) are 
out of scale). 



taken an expansor of 1/v ab in all iterations from the beginning, except that it would 
be impossible from a geometrical point of view, since there is no segment of length 
y/ab (or contractor \^ab< 1) in the generatrix p\. This new limit curve, r v/ ^, is half 
"wrinkled" and half "stretched" , and thus has an intermediate dimension 



Remarks 



dim MF (r 6 ) dim MF (r Vab ) dim A/F (r a 



1. To start at the left (right) endpoint by 1/6 (1/a) is arbitrary: it suffices to take 
any endpoint of any segment with length b (a) to guarantee the polygonal p' k is fitted 
inside p' k+ i after the enlargement, which ensures the inheritance. 

2. Of course, one can expand making another choice, for example: one third of 
the time by 1/b and the remaining two third by 1/a. In this case, for k large, 



Afc ~ (l/6)3 fc (l/a)3 fe = (yf\/a 2 b} , with an expansor c 1 = l/\/a 2 b, c being strictly 
between \fab~ and a. Thus, we have the following proposition. 

Given two contractors a and b, we call p a,b the corresponding generatrix pi of F. 
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Proposition 5.1 For allcGzWL, b<c<a, there exists (up to a translation) a unique 
limit curve T c ^Tmf constructed (not in a unique way) from p a,b . Also 

dim MF (r fe ) < dim MF (r c ) ^ dim MF (r a ) . 



fb\ x 

Proof. Let c be such that b < c< a. The function g(x)=a — is strictly decreasing, 

\aj 

with g(0) = a and g(l) = b. Hence, there is a unique A, < A < 1 such that g(X) = c, 
which satisfies 

b x a l ~ x = c . 

Let {r k } k ^i be a sequence of natural numbers, 0^r k ^ k, such that lim — = A (e.g., 

k~ >+oo k 

r k =[X-k]). Then, we state an expansion process thus: in step k we expand r k times 
by the factor 1/b , and k— r k times by the factor 1/a. Clearly, there is no unique way 
of doing this. The diameter A k in step k will be 

T\ ^ k fl\ k /i\ Xk /1\ (1 ~ A)fc fl\ k 

A fe 



(where a k ~ (3 k means a k j '(3 k — > 1). Let p c k be the corresponding polygonal, thenp£ +1 
inherits pi in a natural way, due to the construction, which guarantees the existence 
of a limit curve T c . The proof of the inequalities of the dimensions is analogous and 
based on the same arguments used in Theorem 4.1 proved in [7]. 

□ 

In the general case of N contractors < a± ^ • • • ^ aw < 1, the proposition is valid, 
taking ai and aj+i instead of b and a. Therefore, if p a ^^-^ a N j s the generatrix of F, we 
have 

Corollary 5.1 For every CjGM, Oi<Ci<a i+ i, there exists a limit curve T Ci &Tmf 
constructed in terms of p ai, --- ,CLN . Also 

dim MF (r a o ^ dim MF (r c ») ^ dim MF (r^+ 1 ) . 

So, every dimensional pair {dim MF (T ai ), dim MF (T ai+1 )} of the discrete MF spectrum 
is made a continuous dimensional interval [dmi/v#F(r ai ), dimM F (T ai+1 )}. Therefore 

Corollary 5.2 For any ceR, a± < c < a^, there exists a limit curve T c e Tmf 
constructed from p ai, - ,aN , and such that 



dim MF (r ai ) ^ dim MF (r c ) < dim MF (r 

We will call the continuous spectrum [dim MF (r ai ), dim A:fi ?(r ajv )], the MF spectrum 
associated with F. From this corollary it is clear, and a remarkable point, that we 
do not need the intermediate contractors; only the minimal a\ and the maximal a at 
are needed to make the discrete MF spectrum to be continuous. Moreover, from 
Proposition 5.1 and the Corollary 5.2, we have 
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Corollary 5.3 For every ai, i^l,N, there exists c, a±<c<aN, and a curve T c , 
stemming from expanding by \ja\ andl/a^ only, such that diniMF(r ai ) = diniMF(r c ). 

Proof. Follows from Proposition 5.1, with b = a\, a = aN and c = cij. 

□ 



6 Relations to Multifractal Spectra 
6.1 Summary of basic notions 

The self similarity concept can be applied to measures. Let Si, ... , Sn be an IFS in 
M. n with ratios dj G (0, 1), 1 ^ i ^ N. Recall that a Borel probability measure \i is 
called a self-similar measure (SSM) if 

N 
i=l 

where pi > 0, and £)i=i Pi = 1. Hutchinson [8] proved that such measures exist and 
are unique, and in this case F = supp(/x), where F is the invariant set for the IFS. 
Besides, it is known that, provided the IFS satisfies the OSC, all reasonable definitions 
of multifractal spectra of \i coincide [9, 10, 11, 6]. They basically are: the 'coarse' 
spectrum, fc, related to box-counting dimension, the 'fine' or singular spectrum, fjj 
related to the Hausdorff one, and the Legendre transform [6]. 
Briefly reviewing, should F be covered by boxes Bj of e-diameter, then for each Bj, 
we define 

_ log fJ.(Bj) 

a \ B i) - — i > 

loge 

and, for a >0, let N a be the number of boxes with a(Bj) «a, then 

logiV a 

fc[pi) = hm , 

£->0 — log e 

if the limits exists. Hausdorff spectrum fjj is rather more related to a local concept: 
if x£F, then 

q(*)=lim l0g/ ^ (£)) , 

and 

where B x {e) is the ball of radius e centered at x. 
For gGlR and e>0, consider 



(5) 
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where the sum is over fi{Bj) > 0, for Bj in an e-grid of MJ 1 . Assuming the limit exists, 
Tfj,(q) is the L 9 -spectrum of fi. 

The functions fc(ot) and r(q) satisfy r(q) = m.i a {qa — fc{o)}. Assuming differen- 
tiability of the functions (which is true for an SSM /i), if for each q the infimum is 
attained at a = a(q), we have q = f' c (a), r(q) =qa — fc(a(q)) and a(q)=T > (q). 
The Legendre transform fi of r is, by definition 



h(a) = inf {qa-r(q)\ 

— oo<q<+oo 



(6) 



In this situation, q and r are related [12] by 



A? 



E ^ T(q) 



(7) 



i=l 



called the partition function due to its formal analogy with the partition function in 
statistical mechanics. 

Thereafter, we will simply write f '■= fc = Ih = II- Also related to the L 9 -spectrum 
are the Renyi dimensions [15]: let {Bj} a partition of M n induced by an e-grid, and 
/i a probability measure supported on F, let pj = l^(Bj). Then for geR, the Renyi 
spectrum of fi is, by definition 



DM 



i fo&lT.jP 

hm 



log£ 

i lo s(s^: 



lim lim 

e^o?^i q — 1 \oge 

l0g(BUp i Pj) 



lim 



lim 



log e 

log(inf J p i ; 
loge 



para q ^ 1 



para q — 1 



para q = +oo 



(8) 



para q 



-oo 



By Eqs. (6) and (5) it is possible to relate the Renyi spectrum to multifractal f(a), 
and so to note that, for q — 0, D = f(a(0)) =dim feox (F), for q— 1, Dx = f (a(l)) = a(l) , 
and for q = +oo and q = — oo, 



a{+oo) 



and 



D. 



a 



-oo 



(9) 



In statistical mechanics, the partitions of a set are always considered of equal size. The 
L q spectrum r(q) is, then, the free energy of the system described by fi as function 
of the inverse temperature q (see for example [13]). This corresponds to the case 
of having all the contractors a,; being equal, which has been widely studied in the 
literature. In this case, the spectrum (a, / (a)) can be calculated in an explicit form 
[5, 14]. 
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6.2 Relation between the MF spectrum and 
multifractal spectra 

Let us consider the partition of different size induced by the IFS yielding a curve 
F G Th, and /i = (pi, . . . ,pn) an SSM on F. So, if F ilt _ tik — S^o- ■ -oSi k (F), then 



\F ilr .., ik \ = a h - ■■a ik = a\ l - ■ ■ a r £ and n(F ily ^ ik ) = p h - ■ -p ik = p\ 1 -- -p 7 ^, for ]T. 
k. We will calculate the f{a) spectrum in terms of the "frequencies" as the /i is 
distributed among the partitions. 

For xEF, let F k (x) the /c-level set F iu _ tik that contains x, so we have 

log/i(F fc (x)) rilogpiH \-r N \ogp N 

a = hm — = hm — = 

fc^oo log|F fc (a;)| k-+oc n log a x H hrjvlogajv , , 

Ai log Pl + • • ■ + X N \ogp N log (p^ 1 • • -p x N N ) 



Ai log ai + ■ ■ • + X N log a w log (a^ 1 • • • a^f) 



a(Ai, . . . , Xn) 



provided the limits of frequencies A» = lim^oo Ti/k exist, 1 ^ % ^ N, where Ti/k are the 
proportions of pi and ai in each k step, so Yli A* = 1. From this expression a known 
result for SSM can also be obtained 

Remark 6.1 ttmm^c^QWxj where 

. log Pi log ^ 

ttmin = mm and a mM = max 



i<i<at log a, i^i<iv log ai 

Proof. Indeed, it can be easily seen that the critical points of the function a(X%, . . . ,Xn] 

log Pi 

subject to V\. Ai = 1, are (0, . . . , i, . . . , 0), whose values of a are -, 1 ^ i ^ N. 

log ai 



□ 



The ubiquitous Stirling formula (as shown below) used in a standard way, yields 

Y,i Ai log Ai 



/(Ai, . . . , X N ) 



Y,iXi log a t 



Let F G Th as above. Let [c/ min , <i max ] be the MF spectrum associated with F . In 
the remainder of this paper, we relate [d min , d max ] to: (a) the Renyi spectrum; (b) 
the range of f(a) values, and (c) the [a m in, QWx] range, choosing, in each case, an 
appropriate self-similar probability measure over F. In the first two cases, we will 
need to "make" probabilities out of ratios. Since Yli a% — L > 1 we will need to 
"contract the contractors", so that the new sum is unity. There are two "natural" 
ways of doing it: J2i a i = L can be written as J2i a i/L = l, and ai/L defined as the 
new Pi, which is what we do in Case (a). The other way is obvious from the very 
symbol dn, since Yli a d { H = 1, so a^ H would be the new pj, that will be Case (b). 

Remark. In case (a) we shrink each ai to ai/L: it is, exactly, as if we zoom out 
generatrix pi (and p&) until we see its length (from far away) to be 1 instead of L > 1. 
In case (b) we shrink the f{a) spectrum, look at it from far away, until its height is 
1, instead of D > 1. 
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6.2.1 Case (a) 

Proposition 6.1 Let pi — a^jL, then 

d m in = -D+oo and <^max = Dq . 

Proof. Indeed, by (9), it suffices to show that a min = d min 

log pi . \og(di/L) . f logL 
amin = nun = mm — = mm < 1 — 



logaj i logaj i I log a 

1 I log(l/ai)J log(l/ai) 
So, a min = dim MF (r ai ) = d min . 

The other equality is true independently of the chosen probabilities pi. Since F is 
self-similar, it is known that dim^(F) = dim.b ox {F), then, it follows from Theorem 4.2 
that d mSLK = dim MF (T aN )=D . 

□ 

6.2.2 Case (b) 

From Sec. 6.1, we can write 

x e F X(x) = (Ai, . . . , X N ) = X : = a 

2^ Xi logai 

for an SSM \i= (pi, . . . ,p^) chosen. If F\ := {x & F : X(x) = A}, then 

F a = \J F x . (12) 

A /a(\)=a 



(b.l) N = 2 

We have weights p\ and and contractors a\ and a binomial measure, Ai = A 
and A 2 = 1 — A. Therefore 

Alogpi + (1-A) logp 2 

« = a(A) = — — — — , 

Alogai + (1-A) loga 2 

so F a = F x , and f (a) = dim H (F a ) = dim H (F x ) . 

For k, the weight of each segment a r 1 a k T r is PiP 2 _r , of which we have the binomial 
coefficient C(k;r), i.e. 

N k (a) 



Stirling's formula 



k\ ~ k k e k V2wk for k — > 00 
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(where means a^jh^ — > 1), allows us to write 

k\ k\ k k ((T\%( r\H; x 



rj r\ (k— r)\ r r (k — r) k ~ r \\kJ V k 
A the frequency of p's, and 1 — A that of (1— p)'s, so 

log N k (a) AlogA+(l-A) log(l-A) 

t(a) = hm ; = — — ; rr- , 

v 1 k^oo -\ g(a r 1 a k r r ) Aloga! + (l-A)loga 2 
i.e. /(a) = /(A), a function of one parameter O^A^l. 

(b.2) N>2 

Now, for the same value of a there are infinite vectors A = (Ai, . . . , Ajy) such that 
a(A) = a. 

For each k, we have C(k;ri, . . . , rjv), the multinomial coefficient, X^i r « = occur- 
rences of weights Pi 1 • • • of a^ 1 ■ ■ ■ , where 

k ) u ■ g _ f"^V t , 

ri,...,rjv/ ri!---rjy! r i 1 '" r jv v \\A;/ V / / 

using Stirling's formula. Therefore 

,/ />» AilogAiH h Atv log Ayv 

Ai iogai + ■ ■ • + AArlogajv 

Since /(a(A)) = dim H (F\), then, /(a) = sup{/(a(A)) : a(A) = a}. 
Thus, we will use Lagrange multipliers, and maximize the function 

h(X) = f(a(X)) subject to a(X) = a . (14) 

Let g be the auxiliary function, and A £ K 

g(X,A) = h(X)-A[a(X)-a } = (15) 

"Ailogpi H h AArlogpjv 



Ai logai H h AArlogaAr 



- «o 



Ai log Ai j h AjylogAjy _ 

Ailogai H h A^loga^y 

So, for l^i^JV 

= — = 

dX z 

= (log Ai + 1) y^Ajlog a,j - log a^Ajlog Xj - A 

j j 

Then 

= A i lo § a i j ( lo § A * + 1 - A lo ^Pi) ~ lo S a i A J log A J ~ A Yl A J log ^ j 



log Pi^Xj log - log (n^2 Xj log pj- 



and dividing by ^ ■ Xj log aj we obtain 

/ j Xj log \j J2 j i°s Pj 

= log Xi + 1 - A log p.; - log Oj J A ^ J - 

Also 

= implies = a(X) — «o • 
Now, for critical A = (Ai,. . ., Ajv) we have, by (16) and (13) 

= log A; + 1 - Alogp, - logOi(/(ao) - Aa ) l<z<JV 
and subtracting from this, the corresponding equation for Ai, we have 

= log(|) - A1 °g(^) - (/W - A «o) , 

from which we obtain 

/„ \ /M-Aao /• \ A 

T _ T / Ul I ( 

A, — A] 



Pi 



and using ^ Aj = 1 we have 



^1 „/(ao)-Aa A 



/(ao)-Aao A 



Pi 



ai P 



I 1 



from which 



/(q )-Aq ^a 
T l Pi 



E /(a )-Aa A 
u i Pi 



and by (17) 



/(a )-Aao A 
y _ Fj_ 

/(ao)-Aao ^A 

i 



for i = l,. . ., A/". Let fi = fi(A) : = f(a ) — Aa , so we can write 

log Aj = log (a? pf ) - log ^ a? Pj A 

calling it'-=^2j af pf (independent of i) for short, we multiply by Aj 

^ log \ = Ai log (a? pf) - \ log(*) , 
we add on i — 1,. . ., N, and then divide by Ai log a; 

£\ Ai log \ _ ^ ^. Ajlog gj A £\ Ai log ft log(*) 



XliAi log ^ J]i^ lo g a i SiAjlogai Z^Ajlogai 
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that is 

f(a ) = 11 + Aa 



E; Aj log o< 
therefore 

log(*) = , 

hence 

Hence, from (18) we obtain 

A, = afpf = a{^- Aa °p} , (20) 
which we replace in (16) and obtain 

E^flogp, £ af W - Aa >f logp, 
a = a (A) = a (A) = -Lr = =4; „ : . . (21) 

i i 

We have then, for AeR, a unique pair (a(A), /(a(A))). If we are given the weights 
Pi,. . .,pn then, once A is fixed, Ej a^p^ = 1 yield -numerically- f2(A), from which 
we obtain a (A) above, and /(a(A)) = fi(A) + Aa(A), or /(a) = fi(A) + Aa. Since A 
is the Lagrange multiplier, we have 

V/i(A) = AVa(A) , 

and, by the very definition of function h (14), we obtain 

Vh(X) = f(a ) Va(A) , (22) 

and then 

A = /'(a ) , 
provided Va(A) t^O. For each 2 = 1,. . ., iV 

<9Ai <9A; yEj Ajloga iy / 

logp 4 Ej Aj log dj - log a< Ej Aj logpj _ log ft - log a* a 



(E, A, log a,) 2 Ej Aj log aj 



then 

Va(A)=0 a = !^^ Vz = l,...,iV, 

log a; 

that, by Remark 6.1, is only true if a (A) is constant. 
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Note. The proof that A is indeed a maximum for h(X) subject to a (A) = «o is 

in the Appendix. Notice that q = —(a(q)) and r(q) = qa(q) — f(a(q)) identify 

da 

fi(A) = fi(/'(a))=n(g) = /(a)-Aa as -r(q), and 



P' 



is now the partition function. 

Should the contractors be all equal, as in the case of the von Koch curve, we would 
have 

Ai = =^— i = l,...,N, 



and the corresponding values of a (A) and /(a (A)). 
If, instead, the probabilities are all equal, we have 



and by (21) 



A, = -4^T • = !,.. -,JV, (23) 



a? 



a? . 1 



V j log— 1 
« = a(fi) = - V = > ( 24 ) 



and by (13) 



a? , / a? 



/(a) = /(a(fi)) 



E^? lE, ^ 
lQ g(E^ 

a n 

Yl ^ l a lo S a 



2° 

~TT~ • (25) 



We will choose this last case: all contractors have the same weight (lit. the same 
importance, the same value), no contractor is more significant than any other. This 
will allow us to give a thermodynamical interpretation to this case (b). So pi = l/N, 
for i = 1, . . . , N, and we want to find, if possible, f2 min and f2 max such that f(a(fl min )) = 
and f{a>(tt max )) = d 

max • 
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Writing Q = D (D Q = dim box (F)=dim H (F)), by (23) we have 



A ?: = ^^- = af° , i = l,...,N, 



and by (25) 



f(a(D )) = D ^ >— = D , 



for the value of a 



a(D c 



n 3 
log — 



£«?°i°. 



according to (24). So, f(a(D )) = d max = dim MF (T aN ), therefore fi max = Aj- 
Clearly, f2 min is the value of Q satisfying f(a(fl)) = d m i n , that is 

O . S / - A . 

a L min o ■ ^mm ? 

— logaj 

a value that cannot be known by analytical means. Yet, we can prove that it does 
exist. Indeed, let < a\ = ■ ■ ■ = a m < a m+ i ^ ■ ■ ■ ^ < 1, i.e. m contractors equal to 
the smallest one, l^m<N. For short, let f (Q) := f (a(Q)) . Then 

Lemma 6.1 

iim m) = lQ f (1/m) . 

a Q 1 
Proof. JiX i (n) : =—^-= then 

22 j " ' z2j (»;•• "<) 

hm xm = { l,m z = 1 '---' m ' 

n^-oo v y | o i = m+l,...,JV . 

Hence, 

t wrvn n /rm J l/mlog(l/m) i = l,...,m, 
n-*-oo 2 = m + 1, . . . , iv . 

Therefore, (25) yields 

Ei Ai(fi) log(A,(fi)) log(l/m) 



lim f(Q) = lim 
n^-oo n-^-oo 2__,i Aj(iiJ logaj logai 
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□ 

Since /(fi) is continuous, and /(-Do) — D = <i max > d min , then, the existence of f2 min 

log(l Tfb) 

(Q m i n <D ), follows from — < rf min . Let us recall that 

logai 

log di i og: j. 

d min = dim MF (r a i) = 1 + = 1 + — ^- . (26) 

log — log — 



Clearly 



1 i„„/ fllH ha m + a m+ iH Va N 

m 



log — i lo S \ 

logm y ai 



< 



logax , 1 , 1 

& log — log — 

log — W r 



log — log — 

a! ax 

Similarly, if < a\ ^ • • • ^ ciN-m' < dN-m'+i = • • ■ = < 1, i-e. m' ratios are equal to 
ajv, then 

Corollary 6.1 

iim m) = lQ f (1/m/) . 

n-»+oo logaiv 

Since Pi—l/N, we have, by Remark 6.1 

. log(l/iV) logiV , 

«min = a(^-oo) = mm — = - — , and 

i^jv loga, log(l/ai) 

log(l/iV) logiV 

«max = a(S2 +0O ) = max 



lsSisJJV logaj log(l/ajv) 
therefore, 

log(l/m) log(l/m') 

/ (ttmin) = — ; and / (a max ) = — . 

logai logajv 

We will compare fl m i n with other values of Q for which the corresponding f(a(Q)) 
can be calculated. 

With fi = 0we have, by (23) 



Ai = =-r = — , 1 = 1,..., iV 

1 j 
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Figure 6: Diagram of an (a (Q), f(a(Q))) spectrum with its significant values (d m i n > 
D\ in this case). The MF spectrum covers the arc connecting dirriMF(r ai ) (blue point) 
to diniMF(r ajv ) (red point). 



Next, for Q = l we have, by (23) 



L 



N 



by (25) and (24) 
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from which 



D := f(a(l)) = 1 V = 1+ " r>! 



We claim: d m i n ^D, i.e. f2 min ^l. 

Since (without loss of generality) a± ^ • • ■ ^ < 1 , we have 1/aj ^ 1/ai for 1 ^ z ^ N, 
hence 

EQj 1 v "\ d% 1 1 — \ Oj 1 

T log — ^ > — log — = log — > — = log — , 

it i 

so, the <i m i n ^-D follows from (26) and (28). 

Note that D\ ^ d m ; n is not always true: the situation depends on the choice of aj. 
Let, for example, a generatrix of F be 




1 

where, N = 4, a\ = = a and a<i = 03 = b. £V cii = L= 2a + 26 and 2a + 6 = 1 . Let us 
write 6 = a p , p^l (1/2 1/3). 

For p = 2, 2a + 6 = 2a + a 2 = l, then a = v^-l « 0.41421, 6 = a 2 = 3-2^- 0.17157. 
Hence 

dmic = 1 + 1£££ = ! + ■°s(^ 2 ) w 1 08983 

log - 2 log - 

a 

whereas by (27) 

— los 4 — los 4 

£>i = /(a(0)) = j ^ = w 1-048585 . 

- log a + - log b - log a 

So, Di<d min . 

But for p= 1.5, 2a + 6 = 2a + a L5 = l (a^0.382, 6^0.236), we have 

d min ^ 1.147 and £>i « 1.152 . 

So, d m \ n 
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6.3 Interpretation of parameter Q, 

The (a, /(a)) spectrum above, with 1/N) q a~ T = 1, — r(g) = f2, f'(oi) = q = 

log(E,«D ( , , , , . , r x W lfil1 logl/iV logl/iV 
(and we have r(o) given by q{T)), fulfills a m ; n = — , a T 



i t. t- \ ■ ■ - — — ■ - ■ \i / o~ ■ " — ~j iv" / / i ~ ~ —-Him i a inax " -. 

log N log a! logaAr 

and Do=max Q f(a>). Now we shrink the spectrum, by shrinking both axes, horizontal 
and vertical, until max Q / sh (a) = l, where / sh is the shrunk spectrum. In case (a) we 
zoomed out the generatrix of F until we "saw" it to be of unit length instead of L^l; 
now we zoom out the spectrum until we "see" its height to be unity, instead of Doll. 



a min log 1/N log 1/N 



As both axes rescale by 1/D , we have now aQ n = — — = — — = , and 

7 A) D log ai log a? 

max = °^ g • We "contracted the contractors" a, to a smaller af° value. Since / 



" ''JV 

fsh 



loga D ° 



and / have exactly the same shape, the slopes of their tangents at any point have 
to coincide: q = q sh . Indeed we can write 



logiV logiV 



where r has been replaced by t/Dq. The partition function above can be thus rewrit- 
ten 

(1/N)" ^ (1/iV)^ 



(°?°) 



t/D 



with o = q sh and r replaced by t/D : this is the partition function for f sh (a). We 
have now Yli = 1 = ET a ?° ■> so we can invert contractors and probabilities [16], 
obtaining a new spectrum: the inverted of the shrunk of the original f(a), which we 

d/*(a) r 

call f*(a). We have now q* = — = , using the permutation q «-> — r for 

da D 

inverse spectra, with t/Dq in place of r. Now — r = Q, so the new derivative is 
Q/D : the inverse temperature for the entropy and internal energy given by f*(ce) 
and a, according to the thermodynamical interpretation. The range Q : 1 — > D is 
here £l/D : 1 — > 1/D : from the information or entropy dimension, up until slope 
1/D (see Fig. 7). ' 

Note. The same thermodynamical formalism applies to Renyi dimensions D q : from 
their definition we have (q — l)D q = r(o), where q is read as inverse temperature for 
the free energy r{q)/q. The abrupt change in the function D q at a certain value of q 
is interpreted as a phase transition -at that q- in many phenomena. 

6.3.1 Case (c) 

Let us now return to our b< y/ab< \/ a?b<a case above. We stressed that, though b 
and a were contractors, \fab was not. It can be written as a^b? and, expanding by 

1 1 ( 1 \ 1 

we would have p' x as expanded by , p' 2 as p-z expanded by 



fab Vab \VabJ ab 
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/*(«*) 



D* = l 




Do 



mm a,. 



a 



Dn 



max (J 



tt' 



Figure 7: Diagram of inverted of shrunk spectrum (a*, /*(«*)) 0/ i/ie one m Fzg. 
i/ere q* can be interpreted as the slopes of tangents to the curve (q* = f*(a*) ). 



p' 4 would be p 4 expanded by [ — ) (in any order: it could be -777 - • • •) etc. If we 

\abj ab 2 a 

restrict ourselves to the even p' k , we have the very same limit curve: T ab , where ab 

would not be a contractor in pi, but in p 2 : we would start from p 2 onwards (i.e. 

skipping the odd pk in the case c _1 = ). With a suitable skipping we can obtain 



2 , 1 

3 



r^ a 6 , which corresponds to as b 
c = a Xa b x \ A a + \ b =l. 
Now, in the general case of a\ 

A 



ab 

and the same is valid for any configuration where 



(Ai 



a^r contractors, let us choose a critical vector 
logl IN 

Xn) corresponding to a certain a = = = , such that 

Ai log aiH \-\ N loga N 

/(A) is precisely /(a), due to the maximization process in Subsec. 6.2.2 (b.2). Since, 

by Eq. (10) 

-log TV logiV 



a 



log 



1 



then, a growing implies log [1/(0^ • • -a^)] decreasing. Notice that a Xl ■ ■ -a^, pro- 
vided that a Xi are rational numbers, is a contractor c that appears for the first time 
in a precise pk and, with a suitable skipping process, starting from p' ko with expansor 



(af 



■ a 



N 



) 1 , generates the curve r a i 1 "' a N N (analogous to Til above). 



For instance, the A critical for a min has the "signature" Q, |, |, |,0), which ensures 
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only o's and no a's (i.e. ai's and no a2's). The signature for afr (square of vab) is 
(I' I' I' I' §)' wn i cn ensures half of o's and half of a's, exactly. The signature for a max 
is (0,0,0,0,1), etc. 

Let us fix a certain critical vector A, to which corresponds a certain (fixed) value of 
a, we are in p k approaching F (recall that any other than A signature (Ai, . . . , Xn)) 
yielding the same value of a, would not produce /(Ai, . . . , Xn) = /(a), but a smaller 
value, corresponding to f{a'), a' y^a. That is why we will work with critical signatures 
of A for each a). 

Fixing a according to (29) implies fixing the length of segments (a^ 1 ■ • ■ a^) k = 
a ? ' ' ' a 7v V ; Z~2i r i = k, since we are in p k (for simplicity, we will refer to the values 
of a in polygonals p k approaching F). These segments of equal length approxi- 
mate, as k grows, an a-Cantordust dense in F. The dim# of this is precisely f(a). 

~\ X\ Xj\j — — 

Going to the corresponding T = T a i "' a N , A = A(a), we find in the corresponding 
nested p' k (with the same suitable skipping) segments of unit length, they are, exactly, 
those segments in p k of length (a Al ■ ■ ■ a^) = a^ 1 • • • a r ^ when expanded by the fc-th 
power of expansor c _1 = [aj 1 ■ ■ -a^) . As a grows, the expansor (a^ 1 • ■ • a^) 1 
decreases... and dimMp(r A ) increases, much as in the first example, when the ex- 
pansor 1/6 = 1/ai = 4 decreases from fc-powers of 4 to fc-powers of 1/a = 2... and 
dimjv/F increases from d min to d max . But then the number of unit segments in the p' k 
of each T A , as k grows, is a precise function of dim M ^(r A ). Let us recall our simple 
example where (ai, a 2 , a 3 , a 4 , a 5 ) = (|>|,|,|>§)- Let us take p k s and p' k J s with k 
even, so c = \fab~ appears as c 2 = ab. If the underlying diameter is unity we have 
16, 8, and 1 segments of length 1/16, 1/8 and 1/4 respectively. Now we expand by 
(l/a) 2 = 2 2 , corresponding to a m i n : we have 16, 8, and 1 segments of length 1/4, 1/2 
and 1 respectively. Now let us expand by the other extreme value a max , by 4 2 : we end 
up with 16, 8 and 1 segments of length 1, 2, and 4 respectively... And, by expanding 
11 (\ 1\ _1 

by — = c _1 = ( — I =8, an intermediate a, we end up with 16, 8, and 1 segments 
ab \2 4J F 6 

of length 1/2, 1 and 2 respectively: the number of unit segments in p' k , as k expands, 
is a give away of diniMF(A). (Remark: should the A not be rational numbers, a limit 
process based on reasoning of Proposition 5.1 would yield analogous results.) 
We stress that we work with critical A only, a (A) yields a contractor and 
an expansor c -1 , as evidenced by (29); c = c(A). To this contractor we have associated 
an a-Cantordust, dense in F, its dim H is, precisely, f(a). To the expansor c~ x (A) we 
have a family p' k converging to a T x e Tmf- T a has unit segments given by expansor c _1 
applied to contractor c in some p k . There is a one-to-one connection thus described 
between with dim AfF : d min — > d max . Two magnitudes grow: a and 

dimjv/F, as two others decrease: the expansor (from 4 fc to 2 fc , in our simple example) 
and the number of unit segments in the corresponding V. The propagation of this 
number of unit segments in p' k as k grows (propagation that can be quantified) is a 
fingerprint of the dim^^ of the corresponding T... But the quantitative law relating 
said propagation of unit segments to the corresponding dim M ^ is, as yet, an open 
question. 
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7 Conclusions 

Contractive processes producing a fractal bounded curve F with Hausdorff dimension 
dim^(F) > 1 can be associated with expansive processes producing locally smooth 
and unbounded curves T with dimMp(r) > 1. To each such F (belonging to an 
ample family of curves) we associate a Mendes France MF dimensional spectrum, 
dimjy./F £ [rfmin, drnax]- Maximal <i max and minimal d m i n of the MF spectrum are identi- 
fied with Hausdorff and division dimensions of F. Such MF dimensional spectrum is 
compared with other multifractal spectra common in the literature, viz the spectrum 
of generalized Renyi dimensions, the thermodynamical formalism (a,f(a)), and a 
one-to-one correspondence between the f(a) dimensions and the dim^p ones, via 
the critical frequency vectors A, which poses an open problem. The MF spectrum 
and the thermodynamical one are compared in terms of universal indices, together 
with their physical interpretation. 



A Appendix: Maximality of A 

To see that the critical frequencies A i; i = l, . . . , N, obtained in (20), are indeed the 
maximum of h(Xi,. . ., Xn), we will use a classical tool from the theory of real valued 
functions of many variables: a determinant called the Bordered Hessian, used for the 
second-derivative test in certain constrained optimization problems ([17]): 

Theorem A.l Let h : U C R N -> R and a : U C R N -> R be of class C 2 . Let XeU, 
a(X)=ao and let S be the level set of a of value ao- Suppose that Va(A)^0 and that 
there exists a real number A such that V/i(A) = AVa(A). Let g(X) — h(X) — Aa(X) be 
the auxiliary function and the bordered Hessian 






da 
<9Ai 


da 
dX 2 


da 
dX N 


da 


d 2 g 


d 2 g 


d 2 g 


<9Ai 


dX\ 


dXidX 2 


dXidXx 


da 


d 2 g 


d 2 g 




dX 2 


dX\dX 2 


dX\ 




da 


d 2 g 


d 2 g 


d 2 g 


dX N 


dXxdX N 


dX 2 dX N 


dX% 



(30) 



evaluated at A. Then, if the determinants of diagonal-submatrices of order ^3 start 
with a subdeterminant \H 3 \ >0, and they alternate their signs (\H^\ <0, \H 5 \ >0,... ; 
etc.), X is a local maximum of h constrained to S. 

Proposition A.l The vector A = (Ai,. . ., Ajv) obtained in (20) is a maximum of 
h(X) (13) constrained to a(X)=ao, where a(X) is the function form Eq. (10). 
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Proof. Clearly, h(X) and a(A) are functions of class C 2 on U Cl*, for 
U = {(Ai, . . . , X N ) ■ X t >0, i = l,...,N] 



and A = (Ai,. . ., Ajv) fulfills XeU. We have also seen that Va(A) 7^0, and for A = f'(a ) 
(22) 

V/i(A) = AVa(A) . (31) 

Partial derivatives of a(X) 



9 m 
dX, a{X) 



d f AjlogPj^ _ log ^ Yij A i lo S a i _ lo S a i Ej Aj \ogpj 
d\ ^EjAjlogaj 



(EjAjloga 



log ^ - log ai 



EjAjlogaj 

EjVogPj 

E 7 - Xj log a* 



log pi - log ^ a(A) 



Then 



9a - _ log p-i - log ai a 
dX i EJ.doga, 



(32) 



O 2 



dXkdX, 



a(A) 



sa; 



Ej V°gaj 



log pi - log ai a (A) 



-loga fc / \ - log ai 9a 

— 2 (log pi - log ai a(A)l + - 



(Ej Ajlogaj 



E, A, log a,) ^ ' 



using (32), we obtain 
9 2 



a(A) 



loga fc da logai 9a 



^^Ai - ^ . A . log a ^ dX, ^ A . log a ^ <9A fe • 

In particular, for k—i 



(33) 



<9A; 



a(A) 



2 log ai 9a 



E.A.loga,) 9Xi 



(34) 



28 



Partial derivatives of h(X) 



_d_ 



h(X) 



d ( Ey \ log Ay \ (log A< + 1) (Ey Ay log a}j - log ai £\ Ay log x i 
dx i \ Ey Ay log ay 

Ey Ay log ay) 



E^jiog^ 
1 



For k^i 

d 2 
d\ k d\ 



h(X) 



Epilog a. 







E 7 - Ay log Ay 

, (log Ai + 1) - log a, 

Ey Ay log ay 

(log Aj + 1) - logaj h(X)j . 



(35) 



dX, 



(Ey Ay log a 
log a k 



(log Ai + 1) - logG^ h(X) 



Ey A y log ca- 



using (35), we obtain 
d 2 



dX k dXi 



h(X) 



(log Ai+1) - log aih(X)) + 



log a k dh 



— log a, dh 



Ey Ay log ay 



dXu 



log a,- 9/i 



Ey Ay log ay) dX i (Ey Ay log ay) 9A * ' 



(36) 



For h — i 
& 2 



dXi 



MA) = ^ 



Ey Ay log ay 
- log Oj 



(log Ai + 1) - log ai a(A) 



2 , (log Ai+1) - log ai h(X) 
Epilog a,) 7 (E, A, log a, 



again, by (35) 



dX] 



h(X) 



1/Ai 



Ey Ay log ay) (E.; A, lo 



2 log ai 9a 
~ dXi 



1 , dh 
A-" l0ga ^ 



(37) 



Second— order partial derivatives of g{X) 

For ky^i, g(X) — h(X)— Aq(A), then, by (33) and (36) 
d 2 g 



d 2 h d 2 a 
dXkdXi dXkdXi dXkdX 
loga fc 



Ej Ay log a 



dh da 
dXi dXi 



log ai 



Ey Ay log a 



dh , 9a 
- A- 



<9A fc 9Aj 
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therefore, from (31) 

For k — i, from (34) and (37) 
d 2 



d 2 9 rr 



d\ k d\i 



(A) = 



(38 



d\ 2 



d2 MA)- - d2 

l/Ai 2 log a; 



MA)-A^a(A) 



(Ej A i lo s (Zj A i lo § %) 



9/i 9a 
<9A,- <9A,- 



then, from (31) 



— (A) = 

5A i ZjAjlogOj 



< 



(39) 



da — d 2 g — 

Next, writing Ai = Tpr-(A) and Bi — -jjttjCA), the bordered Hessian (30) evaluated at 



A, is 








-A, 


-A 2 ■• 


• ~A N 




-A, 


By 


•• 





\Hn+i\ — 


-A 2 





B 2 ■■ 







-A N 





•• 


B N 



(40) 



Notice that, since Va(A) ^ 0, it can be assumed, without loss of generality, that 
Ai^O. Then, for n = 2 








-A, 


-A 2 




\H 3 \ = 


-A, 


Bx 





= B 2 




-A 2 





B 2 





-A 1 
-A, B x 



-A 2 ) 



-A x B L 
-A 2 



-B 2 A\ — B\A\ = B 2 \H 2 \ — B\A\ 



hence, from (39) 



\H,\ > . 



It can be easily seen that, for 3^n^N we have 



n+l\ 



B n \H n \ 



4 )n-l A 2 



B, 



B n -i 



n-l 



B^I + C-l^-^JjBi 



i=l 
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and, by a recurring process 

egn(\H n \) = (-1)" +1 . 
Hence, from Theorem A.l, the vector A maximizes h(X) constrained to a(X) = aQ. 
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